Twin peaks in rf spectra of Fermi gases at unitarity 
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We calculate the radio-frequency spectrum of balanced and imbalanced ultracold Fermi gases in 
the normal phase at unitarity. For the homogeneous case the spectrum of both the majority and 
minority components always has a single peak even in the pseudogap regime. We furthermore show 
how the double-peak structures observed in recent experiments arise due to the inhomogeneity of the 
trapped gas. The main experimental features observed above the critical temperature in the recent 
experiment of Schunck et al. [Science 316, 867, (2007)] are recovered with no fitting parameters. 
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Pairing and superfluidity play a central role in the 
physics of interacting Fermi gases spanning a wide range 
of energies from neutron stars and quark-gluon plasma 
over conventional superconductors to cold atomic gases. 
For weakly-interacting systems the two phenomena ap- 
pear simultaneously, whereas for strongly correlated sys- 
tems pairing sets in at temperatures T* higher than the 
critical temperature for superfluidity Tc. This raises in- 
triguing questions concerning the nature of the paired 
state for a strongly-interacting system. Radio-frequency 
(rf ) experiments on cold atomic gases at unitarity enable 
the study of such fundamental questions P, 01 ■ In partic- 
ular, recent experiments where superfluidity is quenched 
by a large population imbalance examine the strongly 
correlated normal state 0, Q- At low temperatures, the 
shifted single peak spectrum was attributed to the pres- 
ence of pairing without superfluidity. For intermediate 
temperatures, a double-peak structure of the rf spectrum 
was observed. It was speculated whether this double- 
peak comes from the local co-existence of paired and 
unpaired atoms as in the BEG regime. For a balanced 
system at very low temperatures, in-situ measurements 
indicate that no such double-peak structure exists for a 
homogeneous system in the superfluid phase Q]. 

A theoretical description of a strongly-interacting 
Fermi gas is challenging, in particular at nonzero tem- 
peratures In the present paper, we present such an 
analysis based on a diagrammatic approach that includes 
finite temperature pair correlations in the normal phase. 
At the single-particle level, the theory is shown to com- 
pare favorably with Montc-Garlo calculations in the limit 
of a very large population imbalance. We find that for a 
homogeneous system in the normal phase the rf spectrum 
is characterized at all temperatures T and imbalances by 
a single peak, even in the pseudogap regime where strong 
interactions suppress the single-particle density of states 
around the Fermi energy. This holds for both the ma- 
jority and minority components indicating that a sim- 
plified picture of unpaired atoms co-existing with non- 
condensed pairs is too simple at unitarity. Interaction 



effects become more important with decreasing T, and 
the spectra shift to higher energies as compared to the 
bare atomic resonance. For trapped systems, we demon- 
strate that the presence of double peaks in the spectrum 
is due to the inhomogeneity of the sample. In the large 
imbalance case, we recover the main features of the rf 
spectrum observed by the MIT group [1] . 

We consider a gas of fermionic atoms with mass m in a 
potential ^(r) = m(w^a;^ -|-w^2^)/2. Initially, the 
atoms are prepared in a mixture of two hyperfine states 
|1) and 1 2) with an interaction characterized by the s- 
wavc scattering length ai2 . An essentially uniform radio- 
frequency beam couples the state |2) to a third hyperfine 
state 1 3) which initially is empty. We describe the effect 
of the rf beam on the atoms with the operator {Ti — 1) 



rf 



dr 



e-''^*V3(r,t)V'2(r, t)+h.c. 



(1) 



where ■ipi{Y,t) is the field operator for the atoms in state 
f2 is the Rabi frequency describing the coupling of the 
involved hyperfine states to the electromagnetic field, and 
Lo is the rf frequency. The induced transition rate R{uj) 
from state |2) to |3) is within linear response given by 
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where T>{t,y' is the Fourier transform of the re- 
tarded pseudospin flip correlation function —i6{t — 

t')([V'l(r,i)V2(r,t),V|(r',i')'/'3(r',t')])- 

For the two hyperfine states of ^Li used in the MIT 
experiment, a broad Feshbach resonance with |ai2| — > oo 
is located at 834G. Glose to resonance there are strong 
correlations between the atoms in states |1) and |2). Fur- 
thermore, the state |3) in general interacts with the state 
|1) This makes a consistent calculation of the corre- 
lation function D a complicated problem. In the present 
paper, we use a many-body theory which includes the 
two-particle physics leading to the Feshbach resonance. 
This corresponds to treating the interaction between |1) 
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FIG. 1: (a) The atomic propagator for the |2) atoms: thin 
(thick) lines are non-interacting (interacting) Green's func- 
tions, E is the self-energy and F is the T-matrix treating the 
interaction with the |1) atoms in the ladder approximation, 
(b) The correlation function T). 



and 1 2) atoms in the ladder approximation as discussed 
extensively in the literature Q • The diagrammatic struc- 
ture of our calculational scheme is shown in Fig. [TJ We 
neglect for simplicity the interaction between states |1) 
and 1 3). The correlation function T> then decomposes into 
a propagator of |2) atoms interacting with the |1) atoms 
and a free propagator of |3) atoms. Assuming that the 
population of state |3) remains negligible throughout the 
experiment, the transition rate ([2]) for a homogeneous 
system can be written as 



Iml?(tj) = -V 



(2^ 



■^2(fc,6fc-M/(C2fe-0^) (3) 



with V the volume of the system. Here ^ifc ~ k"^ /2m — /i^ 
is the kinetic energy of the \i) atoms relative to their 
chemical potential fjii and /(e) = [exp(/3e) -I- All 
propagators arc calculated within the imaginary time for- 
malism and Wick rotated iw„ —> w + iO+ to obtain the 
retarded correlation function 'D{ijj) at the relevant rf fre- 
quency. We measure frequencies relative to the |2)-|3) 
resonance in the absence of atoms in state |1). The for- 
malism is described in detail in Ref. 0] and here extended 
to the case fJ-i fJ-2- This allows us to model experiments 
with population imbalance S = {Ni — N2)/{Ni + N2), 
where Ni is the number of \i) atoms in the trap and 
Ni > N2. 

We treat the trapping potential in the Thomas- 
Fermi approximation calculating the rf spectrum for the 
trapped system by integrating over the trap the homo- 
geneous signal obtained from Eq. ([3|) with local chemical 
potentials /ii(r) = Hi — V{r). The correlation function 
I?(w) obeys the sum rule 



duj 



ImP(w) = N2 



(4) 



The global chemical potentials /.Ji and ^2 must be ad- 
justed such that they reproduce the correct A^i and 
A^2- The density profiles of atoms |1) and |2) may 
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FIG. 2: (Color online) Rf spectrum of a homogeneous Fermi 
gas with equal populations. Inset: spectral function for a bal- 
anced system at T = 0.3Tf and for the minority component 
of an imbalanced system at T = O.OSTf- (fc^ = y^2mn). 



be obtained from the spectral functions as nj(r) = 
J dk J dLoAi{k,uj)f{uj)/{2TT)'^. Also from this expression 
one may calculate iVi and A^2- For large imbalance, it 
is a good approximation to take the ideal gas values for 
ni(r) and ni since Ni ^ N2. We have checked that 
determining N2 from the interacting density profile and 
from Eq. ([4]) yields consistent results. 

Consider first the rf spectrum of a homogeneous and 
balanced system at unitarity. The system is in the normal 
phase for T > T^, where ksT^ « 0.26Ef from the Thou- 
Icss criterion with the thermodynamic potential calcu- 
lated in the ladder approximation The Fermi en- 
ergy is defined as Ep = (GTT^n)'^^^ /2m, with n = ni = n2- 
We determine the chemical potentials = l-J-2{T) 
sclf-consistently in order to keep the density n fixed. 

In the inset of Fig. [5] we plot the spectral function 
A2{k,Lu) = -2lm[G2{k,uj)] &t T > Tc where the sys- 
tem is in the pseudogap regime: at low momenta A2 
has the characteristic double-peak structure, leading to 
a suppression of the density of states around the Fermi 
energy [1, 0, [l^ . To check the accuracy of our approach, 
we also plot in the inset of Fig. [2] the spectral function 
A2 in the case of a single |2) atom in a sea of |1) atoms 
for very low T. We see that there is a narrow quasi- 
particle peak at w = —Q.QEp, in full agreement with a 
recent variational quantum Monte Carlo calculation [ist . 
This indicates that the ladder approximation used in the 
present paper accounts in the extremely imbalanced case 
for most of the correlation energy at unitarity confirming 
the accuracy of our results [14 1. 

In Fig. [5] we plot the rf spectrum as a function of tem- 
perature for the balanced case. Wc sec that the integra- 
tion over momenta in Eq. ^ has washed out the double- 
peak structure of the spectral function, and the spec- 
trum always consists of a single peak. In the classical 
regime, the |1)-|2) cross section scales as a on \\ with 
At = yj2TT /mksT the thermal dcBroglie wavelength. 



FIG. 3: (Color online) Rf spectrum of the majority and 
minority components (respectively, in the main figure and 
in the inset) of a homogeneous gas for various imbalances 
a = [ni — n2)/{n\ + 712) at T = Q.ZTf and constant ni. 



FIG. 4: (Color online) Rf spectrum of the minority component 
of a trapped Fermi gas with population imbalance 5 — 0.9. 
Each spectrum satisfies the sum rule Eq. (j4| , but here curves 
are normalized such that they reach 1 at their highest peak. 



The scattering rate thus becomes suppressed with in- 
creasing temperature and the system approaches ideal 
gas behavior, exhibiting a narrow spectral line centered 
at the atomic frequency lo = 0. With decreasing tem- 
perature the rf line smoothly shifts to higher energies 
due to increasing pairing interactions between |1) and 
2) atoms, but no double-peaks appear in the spectrum 
of a homogeneous system above Tc, even though the sys- 
tem is in the pseudogap regime with a double-peaked 
spectral function. Similar asymmetric spectra with sup- 
pressed weight a.t u! ~ were observed for the balanced 
system at lower T in the superfluid phase i, m. In 
Fig. [3] we plot the spectrum as a function of imbalance 
for fixed temperature and density ni [l^. The offset of 
the peak for one component is mainly determined by the 
density of the other component. The majority spectrum 
smoothly changes from a broad shifted peak to a narrow 
unshifted peak, in accord with the fact that the majority 
component at high imbalance is an ideal gas. The mi- 
nority spectrum instead remains broad and its offset is 
roughly independent of imbalance. Wc have verified that 
this physical picture holds at all temperatures down to 
T = in the region where the gas is in the normal phase. 
This indicates that the picture of molecules mixed with 
unbound atoms no longer holds at unitarity where the 
pairing is rather a many-body phenomena. By compar- 
ing the width of the spectral peak at unitarity, which is of 
order kp/m, with the molecular energy, of order 
we expect the co-existence picture to become valid on the 
BEC side of the resonance when kpa"^!. 

We now consider the trapped case. In Fig. S] we plot 
the rf spectrum of the minority component for different 
temperatures at imbalance 5 = 0.9 corresponding to the 
experimental conditions of Figs. 2(a-d) in Ref. [3|]. The 
Fermi temperature is defined as ksTp = (6iVi)i/36j, with 
Lj^ = ujx^y^z- Within our theory, the gas does not be- 
come superfluid at this high imbalance. At high T one 
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of the integrand 
0.9 and T = OATf 



FIG. 5: (Color online) Surface plot 
r^ImD(r, oj) of the trapping average at 5 = 
for the minority component. In the inset, we plot n\{r) /n\(Q) 
(dashed), n2(r)/ni(0) (continuous) and r^n2(r)/i?y^n2(0) 
(dot-dashed). 



recovers an ideal gas single peak centered at = 0. The 
spectral line is broadened by pairing correlations, and its 
width is comparable to what is experimentally observed, 
see Fig. 2(a) of Ref. 0]. As the temperature decreases, 
a broader peak emerges from the atomic one; one peak 
is narrow and centered at w = for all T whereas the 
other peak moves to higher frequencies with decreasing 
T, in agreement with the experiments [1, 0|. To anal- 
yse closer the origin of this double-peak structure, we 
plot in Fig. [5] for r = OATp the integrand of the trap- 
ping average r^ImD(r, w), where with T){r,uj) we mean 
the correlation function evaluated with the local chemi- 
cal potentials Hi{r) {uP'r^ = lo'^x'^ + cOy-y'^ + w^z^). One 
sees that for any given radius [corresponding to a ho- 
mogeneous system with densities ni (r) and 712 (r)] , there 
is always a single peak in the rf spectrum in agreement 
with the conclusions obtained above for a homogeneous 
system. For large r the gas is dilute and the spectrum dis- 
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plays locally a sharp peak at the bare atomic frequency. 
This gives rise to the w = peak in the rf spectrum for 
the trapped gas. As r decreases towards the center of the 
cloud, the rf peak shifts continuously towards higher en- 
ergies as the pair correlations between species |1) and |2} 
increase with increasing densities ni(r) and n2{r). The 
integrand contributes most to the total rf signal for radii 
corresponding to the highest value of r^n2(r); as shown 
in the inset of Fig. [5] this corresponds at T = QATp to 
t/Rtf ~ 0.5 with Rtf = kp/mui the Thomas-Fermi ra- 
dius of the majority cloud. For this value of r, the rf peak 
is located &i u) k, 0.5Ep. This is the origin of the second 
peak located at w w 0.5Ep for the trapped system at 
T = OATp. As the gas gets colder, m.axr[r'^n2{r)] shifts 
towards r — and the pair correlations increase leading 
to the pairing peak moving towards higher frequency. By 
varying the imbalance d, we find that the peak position 
is roughly independent of the imbalance, in agreement 
with experiments [l3l and with our findings for a ho- 
mogeneous gas (see Fig. [3]). At very low T the atomic 
peak disappears: the whole cloud is now in the pseudogap 
regime and the spectrum consists of a single broad and 
shifted peak. Note that our theory obtains a remarkable 
agreement with the experimental findings without any 
fitting parameter. It also reproduces reasonably well the 
experimental values for the shift and temperature depen- 
dence of the pairing peak although the peak is predicted 
to emerge at a somewhat lower temperature and with 
a larger shift than what is observed by the MIT group. 
This disagreement could be due to the fact that we have 
neglected interactions between states |1) and |3), which 
are significant at the considered magnetic field. To cal- 
culate the effects of the |1)-|3) interactions requires the 
inclusion of a vertex correction in the correlation function 
V, the so-called Aslamazov-Larkin contribution . On 
general grounds, we expect the inclusion of the |1) — |3) 
interactions to reduce the spectral width and the shift 
of the pairing peak. Sum rule arguments show that this 
shift scales as ^ ^13 for interactions with the same 
high energy behavior, and therefore vanishes in the limit 
0-12 = ai3 [ll, [iOl. The inclusion of the |1)-|3) interac- 
tion in the superfluid phase was recently demonstrated 
to reduce the pairing peak shift [l^ . 

We briefly compare our analysis to other recent works. 



The results in [21|, |2^, |23| are based on a phenomeno- 



logical theory designed to describe pseudogap physics, in 
which an unknown parameter gives rise to a suppression 
in the density of states above Tc- Using this theory, it 
was argued that the double peaks observed in the rf spec- 
trum of a balanced gas are due to the presence of the 



trap 2l| . in agreement with our conclusions. However, in 
[23] it was claimed that a double peak should also exist 
for a homogeneous system, at variance with our findings. 
After the submission of our work, it was argued that the 
appearance of double peaks in the trapped spectrum can 
be traced back to the functional form of the equation of 



state for a fermionic gas [2J]. The approach in [l5| is 
similar to our analysis; they however focus on the T = 
homogeneous superfluid state. 

In conclusion, our theory reproduces the main features 
of the observed rf spectrum for a trapped gas in the nor- 
mal phase, i.e., the unshifted lo — peak which disap- 
pears with lowering T, the emergence of the pairing peak 
leading to a double-peak structure for intermediate T, 
the single pairing peak for low T, and the independence 
of the pairing peak position on the imbalance. We have 
also demonstrated how the recent rf experiments in com- 
bination with a proper microscopic theory provide insight 
into the nature of pairing for strongly-interacting imbal- 
anccd Fermi gases. 

Part of this work was done while two of us (P. M. and 
G. M. B.) stayed at the Institut Henri Poincare-Centre 
Emile Borel. We thank this institution and the IFRAF 
for hospitality and support. Useful discussions with M. 
Zwierlein and especially C. J. Pethick are acknowledged. 
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